Impurity suppression of the critical temperature in the iron-based superconductors 
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We study the impurity suppression of the critical temperature T c of the FeAs superconductors theoretically 
based on the the ±s-wave pairing state of a two band model. The effects of non-magnetic and magnetic impuri- 
ties are studied with the 1 -matrix approximation, which can continuously treat impurity scattering from weak 
to strong coupling limit. We found that both magnetic and non-magnetic impurities suppress T c with a rate that 
is practically indistinguishable from the standard d-wave case despite a possibly large difference of the positive 
and negative s-wave order parameter (OP) magnitudes. This is because the density of states enters together with 
the OP magnitude for the scattering process. 
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Introduction - The alluring prospect of opening a key win- 
dow to understanding the mechanism of high temperature su- 
perconductivity (SC) has attracted fierce research activities in 
the iron based pnictides . The first step towards this goal 
is to establish the pairing symmetry of the FeAs supercon- 
ductors. Many ideas have been put forward to understand the 
seemingly conflicting experimental observations on the FeAs 
materials with regard to the pairing symmetry. Among them, 
particularly appealing is the sign reversing pairing state pro- 
posed by Mazin and coworkers [3]. It is the ground state of a 
two band superconductivity where both pairing order param- 
eters on the two bands have full gaps while acquiring the K 
phase shift between them, which is referred to as the ±s-wave 
pairing state |4|]. It was noticed early on that there is this type 
of solution to a multi-band BCS gap equation a a, and it 
is quite exciting that it seems to be actually realized in the 
pnictide superconductors. A repulsive interband interaction is 
turned to induce pairing by generating the sign reversal be- 
tween the two pairing order parameters. The ±s-wave state 
seems to be able to explain most of the experimental observa- 
tions indicating the full gap behavior [7] as well as a gapless 
behavior 

In this paper, we wish to show that the relative phase of 
K shows up in the impurity suppression of the critical tem- 
perature T c in an interesting way. We employ the T -matrix 
approximation in the weak coupling two band BCS theory. 
In the previous paper, using the same theoretical method, we 
reported that the impurity effects on the ±s-wave state can in- 
troduce an unusual behavior in NMR l/7\ relaxation rate |@]. 
Therefore, it would be interesting to study the effect of impu- 
rity on the T c suppression in this unconventional pairing state. 
The T c suppression by non-magnetic impurity is, as might be 
expected from the sign changing gap nature of the ±s-wave 



state, in between the s-wave and of-wave pairing states. Un- 
expected, however, is that it is indistinguishably close to the 
standard d-wave case despite a large difference of the positive 
and negative OP magnitudes. We also found that magnetic 
impurities are more efficient pair breakers than non-magnetic 
impurities in the ±s-wave pairing state and therefore magnetic 
impurities yield a faster T c suppression rate in the ±s-wave 
pairing state than in the li-wave pairing state although it is a 
marginal difference with realistic parameters. 

Formalism - We study a two band model for the FeAs su- 
perconductors (SC). The details were presented in the refer- 
ence 111 ill . Assuming two SC order parameters, A/, and A e on 
each band, the two coupled gap equations are written as 
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where Vf, e (k ) k l ) is the phenomenological pairing interaction 
originating from the antiferromagnetic (AFM) correlation. 
The above gap equation permits two solutions. When the 
inter-band pairing interaction Vf, e = V e h are repulsive and dom- 
inant over the intra-band interactions, the state where A/, and 
A e have the relative phase of 71, referred to as ±s-wave pairing 
state, is the ground state. The pair susceptibility is given by 
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where N(Q)i,. e are the DOS of the hole and electron bands, 
respectively, and COafm is the cutoff energy of the pairing po- 
tential V(q). 

The impurity effects are included within the T -matrix ap- 
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proximation as 



co„ = co„ +£°(a>„) +E°(co„), 
A h . e = A h . e +Eft(CO„) +E^(a>„), 



(3) 



where co„ = Tn(2n + 1) is the Matsubara frequency, n nnp the 
impurity concentration, and N tot = iVfc(O) +iV e (0) is the total 
DOS. The 1 -matrices T 0,1 are the Pauli matrices x 0,1 compo- 
nents in the Nambu space. The impurity induced self-energies 
are calculated with the 1 -matrix generalized to a two band su- 
perconductivity as jjj, 
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where c — cot So is a convenient measure of scattering 
strength, with c = for the unitary limit and c > 1 for the 
Born limit scattering. (...) denotes the Fermi surface average. 

Because we are interested in determining T c , we take T — > 
T c limit and linearize the gap equation with respect to the order 
parameters. We obtain 
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with N a = N a /Ntot- The pair susceptibility can be written as 
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It is immediately clear that r\ a = 8„ for a single band s-wave 
gap state and there is no renormalization of the pair suscep- 
tibility %a(k) with the impurity scattering. This is just the 
Anderson theorem of T c for the s-wave SC. In our two band 
case, it is more complicated to draw any simple conclusion. 
In particular, the signs of 8/, and 8 e are opposite because of 
the opposite signs of A/, e . 

Before we show the numerical results we can analyze a 
simpler case. The main pairing process in the ±s-wave pair- 
ing state is the inter-band interaction so that we keep only 



Vhe = V e h interactions in the gap Eqs. (1) and (2), and use Eq. 
(fT2l to obtain 
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where X e /f = \/NhNeVheVeh is me effective dimensionless 
coupling constant. This equation can be compared with the 
similarly reduced gap equation without impurities as 
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which yields the standard single band .s-wave result with 
7;° w 1.14co D exp(-l/A. e //). Eq. (O would yield definitely 
smaller T c than r c ° because 8„ is smaller in magnitude than 
r[ a . When both 8„ are set to zeros we obtain another reduced 
gap equation as 
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which is just the case that we would obtain for a double d- 
wave pairing state ifTHl where the anomalous self-energy cor- 
rections (8 a ) are absent because of the sign-changing OP with 
equal sizes. Our case of Eq. < fT~3T > is not straightforward. If 
both 8 a are positive (their magnitudes are always smaller than 
T|o)), the T c reduction would be simply in between the case of 
a s-wave (no suppression) and the case of a d-wave. But in 
the ±s-wave case 8 a will have always opposite signs and as a 
result the T c reduction can be faster or slower than the d-wave 
case of Eq.(15). A simple rule is the following: in the leading 
approximation the reduction rate depends on the sign of the 
quantity (8/, +8 e ). If it is positive, the T c reduction is slower 
than the d-wave case, and if it is negative, the T c reduction is 
faster than the d-wave case. 
We can utilize the relation |AJ 



|A e | = y/KM as T^T C 
found in the minimal two band model in Ref.l 1 1], and obtain 
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From this we can find that (8/, + 8 e ) w (yNh ~ V"e) is al- 
ways positive regardless whether Nh > N e or < N e . There- 
fore, the actual T c reduction should be slower than the d-wave 
case. How much slower will be determined by the magnitude 
of (8/, + b e ) compared to 1 (s-wave limit) and (d-wave limit). 

From the relation (8/, + 8 e ) w (\/~fih — \f^e) 2 we can guess 
that the the T c suppression rate is rather close to the d-wave 
case because the quantity (8/, + 8 e ) is <C 1 unless the differ- 
ence of the DOSs between the bands is unrealistically large. 
In reality, there are more than two bands and also the intra- 
band interactions - which was neglected in the above analysis 
- would make a simple analysis rather difficult. However, a 
practical rule of thumb is that the T c suppression rate by non- 
magnetic impurities in the ±s-wave state should be quite sim- 
ilar to the d-wave case. We will show the numerical results 
obtained by directly solving the gap Eqs. (1) and (2) below. 



3 



both inter- and intra-int. 
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FIG. 1: (Color online) Normalized critical temperature T c /T® vs 
normalized impurity scattering strength r/kgT^ (c = 0). The cal- 
culations are with the realistic bands N^/Ne ~ 2.6 and with the full 
interactions V/,/, , V ee , Vh e and V eh . 



FIG. 2: (Color online) Normalized critical temperature T c /T® vs nor- 
malized impurity scattering strength T/kgT^ (c = 0). The calcula- 
tions are with the realistic bands Ni,/N e ~ 2.6 and with the interband 
interactions V/, e and V e ;, only. 



Magnetic Impurity Case - Let us now turn to the magnetic 
impurity scattering case. For magnetic impurities, if we as- 
sume only an exchange coupling such as S O (where S is 
the momentum of the impurity atom and a is the spin of the 
electrons), we can draw a simple result from the above anal- 
ysis. Because the exchange coupling flips the spin part of the 
singlet wave function 01211 . the result of the magnetic impu- 
rity scattering is to change the sign of A„ in the numerator of 
Eq.(l 1). The final result in the reduced gap equation is to re- 
place 8„ by — 8„ but keeping y\ a the same in Eq.(13). For the 
non-magnetic impurities for the ±s-wave pairing, we had the 
relation 



(l+8 e )(l+5„)>l. 



(18) 



It is equal to 1 for a d-wave pairing state. Eq. ([TBI was the 
very reason why the non-magnetic impurity suppression rate 
in the ±s-wave state is slower than the li-wave case. 

Now, for the magnetic impurity scattering, we have, as dis- 
cussed above 



interband-int. only 
N h /N s = 9 

■ non-magnetic impurity 
— * — D-wave case 
□ magnetic impurity 




1.4 



FIG. 3: (Color online) Normalized critical temperature T c /T® vs nor- 
malized impurity scattering strength T/kgT^ (c = 0). The calcula- 
tions are with the artificial bands N/,/N e = 9.0 and with the interband 
interactions V/ le and V e f, only. 
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It is then immediately clear that the magnetic impurity sup- 
pression rate of T c for ±s-wave pairing should be faster than 
the d-wave case. For realistic parameter values, the suppres- 
sion rates are, as shown in Figs. 1 and 2 from numerical calcu- 
lations, only marginally faster than the d-wave state. In Fig. 
3 we used exaggerated parameter values to demonstrate this 
the point more clearly. This result, however, is only of an 
academic interest because most of magnetic impurity atoms 
would have a much larger potential interaction than the ex- 
change interaction. 

Numerical Results - With the typical band structure of the 
Fe-based pnictides 01311 . we obtained Nh(0)/N e (0) « 2.6 in 



the previous calculations of Ref. lll 111 . With this realistic pa- 
rameter and all interactions included, Fig.l shows the cal- 
culation results of normalized critical temperatures T c /T® vs 
normalized impurity scattering strength T /kgT^ in the unitary 
limit scattering (c = 0). Weaker limit of impurity scattering, 
for example, with c = 1 would just yields twice slower sup- 
pression rate. As can be seen, there are almost no differences 
among all three cases. This is consistent with our analytic 
analysis because (8/, + 8^) w 0.104 <C 1 in this case. Note 
that the normalization of the impurity scattering strength Y by 
r c ° instead of using the gap values A/, e at T = is for conve- 
nience for comparison with future experiments. Also the fact 
that T/k B T c ° - 1 when T c /T c ° -> is a pure coincidence of the 
parameter choice, which is clear in Fig. 3. 
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In Fig. 2, we artificially shut down the intraband interac- 
tions Vhh and V ee ; without the intraband repulsions T c ° itself 
increases by about 40%. But the normalized T c /T® vs the 
impurity scattering strength F/kgT® are indistinguishably the 
same as the case of Fig.l. This result shows that our main 
analytic analysis for the T c suppression with the interband in- 
teraction only will be valid for more complicate multiband 
model in general. 

In Fig. 3, we artificially increase the DOS ratio to N/,/N e = 9 
which is of course unrealistic ratio; this unrealistic parameter 
yields (8/, + 8 e ) w 0.4. The calculations demonstrate that the 
suppression rate of T c indeed follow the trend that we found 
from the analytic estimation. It also demonstrates that in re- 
alistic case the T c suppression of the ±s-wave state by either 
magnetic or non-magnetic impurities should be indistinguish- 
ably close to the case of the standard d-wave SC. 

Conclusions - In summary, we studied the effect of impu- 
rities for the T c suppression on the ±s-wave SC using a gen- 
eralized T -matrix method. The main finding is that despite a 
possibly large difference of the positive and negative s-wave 
OP magnitudes, the T c suppression rate is practically indistin- 
guishable from the standard d-wave case. This is because the 
DOS enters together with the OP for the scattering process. 
As a by-product, we found the subtle difference between the 
magnetic and non-magnetic impurities for the T c suppression, 
which should, however, be a quite small difference in realistic 
case. 
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